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The blazar sequence is reflected in a correlation of the peak luminosity versus peak frequency of 
the synchrotron component of blazars. This correlation has been considered one of the fundamental 
pieces of evidence for the existence of a continuous sequence that includes low-power BL Lacertae 
objects through high-power flat spectrum radio quasars. A correlation between the Compton domi- 
nance, the ratio of the Compton to synchrotron luminosities, and the peak synchrotron frequency is 
another piece of evidence for the existence of the blazar sequence explored by Fossati et al. (1998). 
Since that time, however, it has essentially been ignored. We explore this correlation with a sample 
based on the second LAT AGN catalog (2LAC), and show that is is particularly important, since 
it is independent of redshift. We reproduce the trends in our sample with a simple model that in- 
cludes synchrotron and Compton cooling in the slow- and fast-cooling regimes, and angle-dependent 
effects. 



I. INTRODUCTION 

Active Galactic Nuclei (AGN) with relativistic jets 
pointed along our line of sight are known collectively 
as blazars. This includes those with strong broad 
emission lines, Flat Spectrum Radio Quasars (FS- 
RQs), and those with weak or absent lines, known as 
BL Lacertae objects [BL Lacs;[l|. In general, the spec- 
tral energy distributions (SEDs) of blazars have two 
basic components: a low frequency component, peak- 
ing in the optical through X-rays, from synchrotron 
emission, and a high frequency component, peaking 
in the 7 rays, probably originating from Compton 
scattering of some seed photon source, either inter- 
nal (synchrotron self-Compton or SSC) or external to 
the jet (external Compton or EC). 

Aside from their classifications as FSRQs or BL 
Lacs from optical spectra, subdivided them 
based on their synchrotron peak. They are con- 
sidered high synchrotron-peaked (HSP) blazars if 
their synchrotron peak v s p v k > 10 15 Hz; intermediate 
synchrotron-peaked (ISP) blazars if 10 14 Hz < v^ k < 
10 15 Hz; and low synchrotron-peaked (LSP) blazars if 
v°l < 10 14 Hz. Almost all FSRQs are LSP blazars. 

[3] combined several blazar surveys and noticed an 
anti-correlation between the luminosity at the syn- 
chrotron peak, L*|, and the frequency of this peak, 
Vp k . They also noticed anti-correlations between the 
5 GHz luminosity L5 GHz and z/jf, the 7-ray luminos- 
ity and i/p k , and the 7-ray dominance (the ratio of the 
EGRET 7-ray luminosity to the synchrotron peak lu- 
minosity) and z/fc. [4j provided a physical explanation 
for these correlations [see also [|| . If the seed photon 
source for external Compton scattering is the broad- 
line region (BLR) , and the BLR strength is correlated 
with the power injected into electrons in the jet, one 
would expect that more luminous jets have stronger 
broad emission lines and greater Compton cooling, 



and thus a lower 
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As the power injected in elec- 



thcre are fewer seed photons for Compton scatter- 
ing, and consequently the peak synchrotron frequency 
moves to higher frequencies. This is also reflected in 
the lower luminosity of the Compton-scattered com- 
ponent relative to the synchrotron component as v ■ 
moves to higher frequencies. 
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anti-correlation has been questioned. 
Using blazars from two large surveys, @ did not find 
any anti-correlation between v^f. and radio, BLR, or 
jet power. This work, however, has been criticized for 
its relatively poor SED characterization Q. Also, the 
lack of sources in the upper right region of the L^ v k - 
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plot could be the result of a selection effect 
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lid? ]. Since a large fraction of BL Lac objects have 
entirely featureless optical spectra, their redshifts, z, 
and hence luminosities, are impossible to determine. 
These could be extremely bright, distant BL Lacs that 
would fill in the upper right region. [Tl| found no 
correlation between the frequency and luminosity of 
the synchrotron peaks for objects in the Metsahovi 
Radio Observatory BL Lacertae sample. [l2| did find 
an anti-correlation, using sources found in the LAT 
bright AGN sample 0, [l3|. In recent works, an "L" 
shape in the L 
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trons is reduced, the broad line luminosity decreases, 



plot seems to have emerged, as 
lower luminosity and low-peaked sources have been 
detected with more sensitive instruments [13, HH ■ [H| 
found more of a "V" shape, although their plot did 
not include FSRQs; if high luminosity and low-peaked 
FSRQs were added, it might appear as more of an 
"L" . With the advent of the Fermi Telescope era, it is 
now possible to also characterize the Compton peak 
frequency v^ k and luminosity L pk . As we show, the 
Compton dominance Ac = Lp k /L s ^ k is an important 
parameter for charcterizing this sequence. 



II. THE 2LAC BLAZAR SEQUENCE 

The 2LAC [16[ allows for the characterization of 
the high energy component for a greater number of 
blazars than previously possible. Here we look at the 
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blazar sequence among the 2LAC clean sample, which 
includes 885 total sources, with 395 BL Lacs, 310 FS- 
RQs, and 156 sources of unknown type. 

[2] fit the broadband SEDs of the blazars in the 3- 
month LAT bright AGN sample [LBAS;[l3[ with third 
degree polynomials to determine the peak synchrotron 
frequency, z/^. found empirical relations for find- 
ing the peak frequency of the synchrotron component 
from the slope between the 5 GHz and 5000 A flux 
(a r „), and between the 5000 A and 1 keV flux (a ox ). 
|16| used these empirical relations and their low energy 
data to determine v sy k for the 2LAC sample. This was 
then used to classify the SEDs of the blazars as LSP, 
ISP, or HSP. also provided an empirical formula for 
determining the flux at the synchrotron peak from the 
the 5 GHz flux density and v s 3.. We used this relation, 



along with the vV 1 



values from the 2LAC |1( 
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the luminosity distance and create a plot of the peak 
synchrotron luminosity, L sy k versus u*%. This is plot- 
ted in Fig. [TJ These are the sources in the 2LAC clean 
sample with measured redshifts and enough SED mea- 
surements to determine v^, which amounts to 352 
sources including 145 BL Lacs, 195 FSRQs, and 12 
AGN of unknown optical spectral type (AGUs; i.e., 
unknown whether they are FSRQs or BL Lacs). Note 
that i/pf. is corrected for redshift, and is in the frame 
of the source. We include the photometric redshifts 
which have been determined from [l7j . 

For this diagram, we have computed the Spear- 
man (p) and Kendall (r) rank correlation coefficients 
and the probability of no correlation (PNC) calculated 
from each coefficient. The results can be found in Ta- 
ble HI The results from p and r are similar in all cases. 
The PNC is very small for the BL Lacs and FSRQs 
separately, and even lower for all the sources com- 
bined, where the probability is essentially zero that 
there is not a correlation. Note however, that, sources 
with unknown z are not included. This could explain 
the anti-correlations for the whole sample and for the 
BL Lacs in general fl8| . since they could fill in the 
upper right part of this diagram, as mentioned in Sec- 
tion |H However, the objects with unknown z are al- 
most certainly BL Lacs, so it would not explain any 
possible anti-correlation amoung the FSRQs alone, al- 
though no significant one was found. 

[H found a more significant anti-correlation between 
the 5 GHz luminosity, L5 GHz, and z/|. We plot this 
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FIG. 1: Peak synchrotron luminosity versus peak syn- 
chrotron frequency for objects in the 2LAC clean sample. 
Filled circles represent FSRQs, empty circles represent BL 
Lacs, and filled squares represent objects which do not 
have an unambiguous classification. Dashed lines indicate 
the boundary between HSPs and ISPs and between ISPs 
and LSPs. 




FIG. 2: Radio luminosity at 5 GHz versus peak syn- 
chrotron frequency for the 2LAC clean sample. Symbols 
are the same as in Fig. [TJ 



determined in part based on the 5 GHz flux, this is 
probably a result of this dependence. As pointed our 
by [H|, if v s y. increases, but the synchrotron bump 
remains unchanged in other aspects, the radio flux 
(or luminosity) will naturally decrease. This can ex- 
plain this anti-correlation. However, what is unclear 
is whether a v sy k derived from the radio flux should 
be interpreted as a cooling break since these two 
should be from different regions and possibly indepen- 



for our sample in Fig. [2] Like [3j, we find the anti- dent. Determination of v sy k independent of low radio 



correlation with L§ qhz versus v sy k to be much more 
visually apparent than the one between L s \ and v^f,. 
This is also reflected in the decreased PNC for this di- 
agram (Table HI for both FSRQs and BL Lacs alone. 
If the anti-correlation is explained by the increasing 
cooling at higher luminosities [3] , the correlation with 
Lpl should be more significant, since the emission at 5 
GHz is thought to be from a different region of the jet 
than the emission at the peak. However, since v sy k was 



pk 

frequency emission should be a good, although tech- 
nically challenging, way to test this. 

[2| also fit the high energy components of LBAS 
blazars with a third degree polynomial to determine 
the peak of the 7-ray component (presumably from 
Compton scattering). They found an empirical rela- 
tion between v^ k and the LAT 7-ray spectral index, 
r 7 . We use this relation to determine the peak of the 
Compton component. Approximately 10% of the 352 
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TABLE I: Statistics of correlations involving v s 



Sample 


P 


PNC(p) 


r 


PNC(r) 


t sy sy 

L p l versus vj. 


BL Lacs 


-0.19 


0.019 


-0.12 


0.035 


FSRQs 


-0.12 


0.073 


-0.088 


0.066 


All sources with known z 


-0.54 


4.2 x 10" 28 


-0.35 


0.00 


Lsghz versus i/Jjj! 


BL Lacs 


-0.64 


5.4 x 10~ 18 


-0.45 


0.00 


FSRQs 


-0.36 


1.82 x 10~ 7 


-0.25 


1.79 x 10" 7 


All sources with known z 


-0.79 


0.00 


-0.58 


0.00 


Ac versus v 3y k 


BL Lacs 


-0.30 


2.7 x 10~ 4 


-0.21 


1.3 x 10 -4 


FSRQs 


8.9 x 10" 3 


0.90 


6.3 x 10" 3 


0.89 


All sources with known z 


-0.66 


9.8 x 10~ 45 


-0.45 


0.00 


All sources 1 


-0.66 


0.00 


-0.46 


0.00 


All sources 2 


-0.66 


0.00 


-0.46 


0.00 


All sources 3 


-0.64 


0.00 


-0.44 


0.00 



sources are also found in the 58-month Swift Burst 
Alert Telescope (BAT) catalog. For these sources, 
we extrapolated their BAT and LAT power-laws and 
found where they intersected. If they intersected 
within the range 195 keV to 100 MeV we used this lo- 
cation as fp k . This allows an improved estimation of 

v pk over the empirical relation, particularly for those 
very soft sources, for which this empirical relation is 
untested. For other sources, the peak was determined 
from the LAT spectral index and the empirical rela- 
tion from Q . Once the location of i/jj. is known, either 
by using the empirical relation or from the BAT-LAT 
intersection, the luminosity at the peak, L^ k , can be 
estimated by extrapolating the LAT spectral index. 

Combining their results with EGRET data, [U 
made a plot of 7-ray dominance versus v^. Using 
LAT data from the 2LAC, as described above, we 
make a similar plot, although we use L^ k instead of 
simply the 7-ray luminosity. Our results are in Fig. [3l 
Also note that Ac is independent of redshift. Thus we 
can plot in Fig. [3] an additional 174 sources from the 
2LAC clean sample that have well-determined syn- 
chrotron bumps but do not have known redshifts. For 
these sources, the plotted v s ^ k is a lower limit, since the 
redshifts are not known. However, it will be larger by 
only a factor (1 + z), i.e., a factor of a few. 

We have also computed the correlation coefficients p 
and t for Ac versus v^ k , and the results can be found 
in TableQ] There is no evidence for a correlation for the 
FSRQs alone, although the probability that there is no 
correlation for the BL Lacs alone is very small. For the 
combined sample of all sources with known z, there is 
essentially zero chance that there is no correlation, 
similar to the L s ^ k versus v s p v k correlation. We also 
computed the coefficients for all sources, including the 



ones with unknown z, computing their v s pk assuming 
z = ("all sources 1" in Table IJ). We find essentially 
no chance that the addition of sources with unknown 
z will explain the correlation when these sources are 
included. The term u^ k will vary by a factor of a few 
due to redshift, so we also calculated the coefficients 
assuming all these sources with unknown z are at z — 
0.35, the average of the BL Lacs with known z ("all 
sources 2"); and assuming these sources are at z = 4 
("all sources 3"), which is higher than the maximum 
redshift of the entire sample (which is z — 3.1). In 
each case, there is essentially a 100% chance that an 
anti-correlation exists. Although the objects without 
known redshifts could explain the correlation between 
Lpl an d v pii it does not appear they can explain the 
correlation between Ac and v s pk . This aspect of the 
blazar sequence seems secure. 

The two sources with unknown redshifts in 
the upper right quadrant of this diagram are 
2FGL J0059.2-0151 (1RXS 005916.3-015030) and 
2FGL J0912.5+2758 (1RXS J091211.9+27595) with 
LAT spectral indices of F 7 = 1.15 ± 0.36 and T 7 = 
1.20 ± 0.37, respectively. These are the two hard- 
est sources in the 2LAC, and the sources with the 
largest error bars on their spectral index; only one 
source in the 2LAC clean sample is fainter than these 
sources (2FGL J1023.6+2959). They are clearly out- 
liers. Propagating the error on their spectral indices, 
one finds that they have Compton dominances of 
log 10 (A c ) = 1.44 ± 1.82 and log 10 (A c ) = 1.48 ± 1.98, 
respectively; they have Ac consistent with unity 
within their error bars, and so are consistent with the 
"L" shape seen in Fig. [3] 
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FIG. 3: Compton dominance (i.e., L^ k /L s ^ k ) versus peak 
synchrotron frequency. Symbols are the same as in Figures 
[T] and [2] additionally, rightward-pointing triangles repre- 
sent BL Lacs with unknown redshifts, for which i/^ k is a 
lower limit. 



III. THEORETICAL BLAZAR SEQUENCE 

We describe a simple model for blazar jet emis- 
sion. This model is similar to the one presented by 
0. We assume the relativistic jet is dominated by 
emission from a single zone which is spherical with 
radius R b in its comoving frame, and moving with 
high relativistic speed /3c giving it a Lorentz factor 
T= (l-py 1 / 2 . The jet makes an angle to the line of 
sight 9 — arccos(/i) so that the Doppler factor is 5^, = 
[r(l — /3/j,)] . Electrons are injected with a power- 
law distribution given by Q e {l) = Qal q H{l] 71,72). 
The hard X-ray spectra in some blazars may indicate 
very hard electron spectra at lower energies [20]. Al- 
though blazars can be quite variable on timescales 
as short as hours [e.g., [2l[ or even minutes [22j] . we 
will assume their average or quiescent emission can 
be described by a steady state solution to the electron 
continuity equation, where continuous injection is bal- 
anced by cooling and escape. We assume an energy- 
independent escape timescale given by t esc — r/R' b /c 
where R' b is the comoving size of the blob and 77 is 
a constant > 1. In this case, where 71 < 7c (the 
slow- cooling regime) the electron distribution can be 
approximated as 



O^esc 



7i < 7 < 7c 



v ilhc) q 1 7c < 7 < 72 
If 7 C < 71, i.e., the fast-cooling regime, 



JV e ( 7 ) Pa Q t e 



(7/71 )~ 
(7/71) 



-9-1 



7c < 7 < 7i 
7i < 7 < 72 



(1) 



(2) 



Here we assume the electrons are cooled by syn- 
chrotron emission and Thomson scattering, so that 



3m P c 2 



4ccr T (u' B + u' + T 2 u ext )t e 



(3) 



is the cooling electron Lorentz factor, where primes 
denote quantities in the comoving frame of the blob. 
Here u' B = B 2 /(8ir) is the magnetic field energy den- 



sity, u 



B 
I 

sy.tot 



is the total synchrotron energy density, 



and the u ex t is the external energy density, assumed to 
be isotropic in the proper frame of the AGN. The ex- 
act nature of the external radiation field is not known, 
and may not even be the same for all blazars. For a 
given set of parameters, the nonlinear nature of u' sy 
means that there is not a simple closed form solution 
for N e (j), and so we solve for N e (j) numerically. 

Once -/V e ( 7 ) is known, we use the standard formu- 
lae in the ^-approximation [e.g., [H, [H, [25[ to cal- 
culate the synchrotron, SSC, and EC luminosity and 
energy at their peaks. This assumes the Compton- 
scattering takes place only in the Thomson regime. 
We have compared our results with those using the 
full synchrotron emissivity and Compton cross sec- 
tion, and found good agreement. The Compton 
dominance (Ac) is given by the ratio of the peak 
Compton-scattcrcd component to the peak of the syn- 
chrotron component, Ac = max[ 



EC 



A L P k 



T SSC] IT 

L pk \I L 



pk 



max^QUeat, u' sypk ]/u' B , where we have ignored a 
bolometric correction term ~ 1. Note that SSC emis- 
sion has the same beaming pattern as synchrotron, 
and thus /L s ^ k does not depend on the viewing 
angle; however, L® k does not, and thus L^P j ' L S JL is 
dependent of the viewing angle through <5d [2y, |27| . 

We found we could adequately reproduce the 
sources in Figures Q] and |3] without changing the power 
injected in electrons. Instead, we only vary B and 
u ex t, assuming B oc u\ x \, and the angle to the line of 
sight, 9. Further details can be found in 28]. 

Using the model described above, we reproduce the 
Lpk versus v 8 ^, relation as seen in the curves in Fig. 21 
At < 10 13 Hz (the exact value depends on 8), j c will 
become less than 71 , leading to a sharp inflection. By 
changing the viewing angle, the synchrotron luminos- 
ity decreases dramatically, so that almost all of the BL 
Lacs can be reproduced by this model. Note that this 
model predicts that there will be sources found with 
Lpk ~ 10 47 erg s _1 and v** > 10 15 Hz, which is rela- 
tively unpopulated. This region could be filled in by 
bright BL Lacs at high z, which for which we cannot 
determine a rcdshift due to the nonthermal emission 
washing out their weak emission lines [181 ] . Indeed, 
this region is beginning to be filled in by constraining 
the redshifts of several BL Lacs [29| . 

The model curves for Ac versus v s pk are plotted in 
Fig. [5] Again, the model seems to reproduce the 2LAC 
data well. There is a distinct difference above and 
below the transition between the Compton compo- 
nent being dominated by SSC and being dominated 
by EC, at around 1$ w 2 x 10 15 Hz for the 9 = 
curve. Note that Ac is not dependent on 9 for SSC, 
although the curve here does shift to lower v s pk due to 
its dependence on 9. The EC part of the curves, how- 
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FIG. 4: Same as Figure[T] with curves showing our model 
plotted at various angles, shown in the legend. Along the 
curves only B and u ex t are varyed. 




FIG. 5: Same as Figure [3] with curves showing our model 
plotted at various angles, shown in the legend. Along the 
curves only B and u ex t are varyed. 



ever, are strongly dependent on 6 through (5d- As with 
Fig. HI a sharp inflection is seen at lower frequencies 
( v pk ~ 10 13 H z f° r ^ — 0) due to the transition of the 
peak from being associated with 7 C to 71 . Also as with 
the previous figure, the BL Lacs are well-reproduced, 
while the FSRQs are not, particularly those with high 
Ac- A diversity of 71 would clearly allow a widen- 
ing of this branch, and allow the model to reproduce 
all of the FSRQs. Alternatively, if the jets of blazars 
are made up of multiple emitting components with 
different Lorentz factors, so that at different angles a 
different jet component would dominate the SED, this 



could also broaden the sequence and reproduce more 
objects [HJ ■ 



IV. DISCUSSION 

For blazars, it has been suggested that the peak 
of the electron distribution {^ p k) is independent of 
broad emission line luminosities, B, and other prop- 
erties [lH . If this were the case, and if the broad lines 
are the source for external scattering, sources would 
be found with high Ac {Ac ^ a few) and and high 



v 



pk ( v 7k ~ 1015 Hz )- The redshift is not relevant, 
since both of these quantities are essentially redshift 
independent. However, no such sources are found in 
our sample, including BL Lacs without known red- 
shift (Fig. [5]). If 7 rays are produced by external 
Compton processes, as seems likely [3(|, there does in- 
deed seem to be a relation between u ex t and 7 P fc . We 
have found a simple model that explains this, where 
Ipk = max[7 c , 71], where 7 C is dependent on both u ext 
and B (Equation [3]). This is quite similar to pre- 
vious work [e.g., y, l5|], the main different being the 
dependence of B and u ex t- This simple model pre- 
dicts sources will be found with L S A > 10 47 erg s _1 
and v s p y k > 10 15 Hz (Fig. |J), when the redshifts of 
more BL Lacs are determined [unlike 0, @- These 
sources should be brighter because they are more 
highly aligned with our line of sight, rather than in- 
trinsic brightness. [3l| also predict that sources will 
be found in this region, although in their case these 
sources will be "blue FSRQs" , with the primary emit- 
ting region ourside the BLR, assuming in other FS- 
RQs, broad line photons are the external radiation 
source for Compton scattering. 
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